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Abstract: We study the potential to observe CP- violating effects in SUSY cascade 
decay chains at the LHC. We consider squark and gluino production followed by 
subsequent decays into neutralinos with a three-body leptonic decay in the final step. 
Asymmetries composed by triple products of momenta of the final state particles are 
sensitive to CP-violating effects. Due to large boosts these asymmetries can be 
difficult to observe at a hadron collider. We show that using all available kinematic 
information one can reconstruct the decay chains on an event-by-event basis even in 
the case of 3-body decays, neutrinos and LSPs in the final state. We also discuss 
the most important experimental effects like major backgrounds and momentum 
smearing due to finite detector resolution. We show that with 300 fb~^ of collected 
data, CP violation may be discovered at the LHC for a wide range of the phase of 
the bino mass parameter Mi. 
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1. Introduction 



Supersymmetry (SUSY) |I|, ^, ^ is one of the best motivated extensions to the 
Standard Model (SM). According to the latest fits [|, |H] to the precision observables 
a rather hght scale is indicated and is in reach of the Large Hadron Colhder (LHC). 
If SUSY happens to be reahsed in nature, the task will be to pin down the specific 
characteristics of the underlying model which will require many detailed studies. One 
of the interesting issues in this context is CP violation. While the observed amount 
of CP violation in the K and B sectors can be accommodated within the SM, another 
piece of evidence, the baryon asjnnmetry of the universe, requires a new source of 
CP violation 0, 

The Minimal Supersymmetric Standard Model (MSSM) fig, |TT], g has about 
one hundred free parameters and a large number of these may have non-zero CP- 
violating phases, see e.g. [0. Many of the phases are unphysical in the sense that 
they can be rotated away by a redefinition of the fields but not all can be sys- 
tematically removed in this way. For example in the sectors of the charginos and 
neutralinos, the supersymmetric partners of Higgs bosons and SU{2) x f/(l) gauge 
bosons, we can have up to three new phases. These include the phases of the U{1) 
gaugino mass parameter Mi, SU{2) gaugino mass parameter M2 and the Higgsino 
mass parameter /x. One of them can be rotated away and therefore we choose the 
convention where M2 is kept real. Hence we have 

Mi = |Mi|e*'^S ^i=\^i\e'^^. (1.1) 

The other possible sources of complex parameters in the MSSM Lagrangian are 
trilinear couplings Af but they are not relevant to this study [|1^ . 



Certain combinations of the CP-violating phases are constrained by the experi- 



mental upper bounds on various electric dipole moments (EDMs), see e.g. [|^. Ig- 
noring possible cancellations, the most severely constrained phase is that of /i which 
contributes to the EDMs at the one- loop level and has a tan/3 enhancement. In gen- 
eral for (9(100) GeV masses, |0^| < O.OItt and we therefore set 0^ = throughout 
our study. The phase of Mi has weaker constraints. Large phases can be accom- 
modated in models with heavy first and second generation squarks and sleptons 
(>TeV) or if accidental cancellations occur [^, 0, |19|, Here we study 



the complete range of CP phases in order to see the general dependencies exhibited 
by our observables and the luminosity required to observe these within the LHC 
environment. Nevertheless, we want to stress that in the chosen scenario experimen- 
tal bounds from EDMs can be evaded by arranging cancellations between various 



supersymmetric contributions for any value of 0i |2^, [T^ . 

CP-odd observables are the unambiguous way of discovering hints of complex 
parameters in the model. One class of these observables include rate asymmetries of 
cross sections, branching ratios and angular distributions. The other possibility are 
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the observables that can be defined using triple products of momenta and/or spins 



of particles, see [24, OR 26 1 for a recent review. In any case studies of CP- violating 



phases will be challenging at the LHC [^, [T^, |2^ and a precise determination of 



them is expected to be only possible at a future e^e~ linear collider. 

An interesting example of using triple product correlations is studying the 3- 
body decay of a neutralino X2 HI) S' EH Since the asymmetry is maximal in 



the rest frame of the decaying neutralino, the central idea of ||2^, |3^ was the full kine- 
matic reconstruction of the production and decay process. Hence the CP-sensitive 



observables are calculated |2^, ^ in the rest frame of the decaying neutralino at an 
e^e~ collider. This can be done due to the clean experimental environment and the 
well-known initial state at such a machine. As this is not the case at the LHC, these 
observables are much more challenging due to the proton structure and further ex- 
perimental uncertainties, hence making the observation of CP-violating effects very 
difficult. However, in the present analysis we show how a similar approach can be 
applicable at hadron colliders. By applying full event reconstruction, the observation 
of CP-odd asymmetries becomes feasible. 

In many scenarios the highest cross section at the LHC is predicted for the 
associated production of gluinos and squarks of the first two generations. If squarks 
can decay to a neutralino X2 process 

qi ^xl + q (1-2) 

one can expect an abundance of neutralinos. Large statistics will help to overcome 
the limitations of a similar analysis of the stop carried out in [l^. The production 



of the neutralino X2 is followed by the subsequent three-body decay 

Xl - xTr. (L3) 

The 3-body decay is discussed since sizeable contributions to CP violation are ex- 
pected from interference diagrams and it offers the opportunity to use CP-sensitive 
observables based on triple product correlations. The process is sensitive to any 
phase that enters the neutralino sector (0i and 0^), but we set 0^ to zero as dis- 
cussed above and concentrate on the effect of 0i. To extract information about the 
phase of Mi we analyse the following triple product of momenta from the decay 
products of the g^, 

T = P,-imy<Pi-). (1.4) 

It originates from the covariant product present in the amplitude and can be ex- 
panded in terms of explicit momenta, 

^/.upaP'^PbPcPd ^ EaPt-ilfcy<Pd)±--- (1-5) 

When CP phases appear a non-zero asymmetry can be constructed from final state 
momenta using the triple product in Eq. ( |1.4| ) as explained below. 
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The triple product T changes sign under the naive^ time reversal T^v and con- 
sequently it is a TTv-odd observable. Neglecting possible absorptive phases due to 
loops, Breit-Wigner propagators and final state interactions, that are expected to be 
small in the leptonic neutralino decay, CPTjv invariance is equivalent to CPT invari- 
ance |^ and hence a measurement of a T^r-odd observable implies CP violation. In 
addition as these are CP-odd quantities any measurement of the phase made by this 
method not only reveals the magnitude of the phase but also the sign. It is important 
to note that the triple product correlations are a tree-level effect and therefore large 



asymmetries can be expected ||32|, p6| . 

To be able to successfully observe CP violation we need to isolate squarks from 
antisquarks as the CP-conjugated process will have an asymmetry of the opposite 
sign. Normally this is done via charge identification but this will be impossible using 
the above decay chain with the light quark in the final state. However, for the first 
two generations, squark production typically dominates over anti-squark production 
due to the valence quarks present in the proton. Consequently our sample will have 
a majority of squarks (around 80% in the analysed scenario) and the asymmetry will 
not be washed out. The proportion depends on the masses of squarks but this does 
not affect our analysis significantly. 

One of the major problems when trying to measure CP-sensitive observables 
with triple product correlations at the LHC are the large undetermined boosts given 
to produced particles. CP asymmetries have the largest value in the rest frame 
of the intermediate particle and the effect of the boost produces a severe dilution 
factor. To remove this dilution factor we have to boost all the particles in the triple 
product {q,i^,i~) into the rest frame of the intermediate particle in the decay (xg 
in our case). In SUSY decay chains this is complicated due to the presence of two 
LSPs and possible neutrinos that escape detection. Once the missing momenta are 
reconstructed it is then trivial to find the momentum of the particle that we are 
interested in and boost back into its rest frame. In principle this should return the 
CP asymmetry to the magnitude one would expect if the initial particle was produced 
at rest. 

So far in the literature the momentum reconstruction was used to determine 
the unknown masses of SUSY particles in the case when the decay chain proceeded 



via two-body decays [0, In this analysis we change the approach and 

assume that the masses are already known with a reasonable precision from other 
measurements at the LHC. The new complications are three-body decays of super- 
symmetric particles (charginos and neutralinos) and possible missing neutrinos which 
means that we have less on-shell constraints. However, as we show in this paper, 
even in such a case momentum reconstruction is still possible after the inclusion of 



^Nai've time reversal T^v reverses the momenta and spins of all particles without interchanging 
the initial and final states, see e.g. [Q. Recall that under the true time reversal also the initial and 
final states are interchanged. 
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experimental uncertainties. Having all the momenta reconstructed one can perform 
a more detailed analysis of particle properties which could open a window for a more 
accurate determination of the underlying model at the LHC 

Finally, in order to assess the feasibility of completing this study in the demand- 
ing experimental environment of the LHC we perform a simulation of our process 
using the Monte Carlo event generator Herwig++ pi] , |4^ . We apply simple selection 
cuts and smear final state particle momenta to account for finite detector resolution. 
Moreover, we make an analysis of the impact of the most important background pro- 
cesses relevant for our study. All of the analysis was done assuming a centre-of-mass 
energy of 14 TeV but it would not change markedly if we moved to 10 TeV except 
for an increase in the luminosity we require to see a statistically significant result. 

The paper is organised as follows. We begin by describing the process under 
consideration in Section ^ and discuss the triple product correlations in more detail. 
In Section ^we present the method we use to reconstruct the momenta of the invisible 
particles in our decay chain. Section || discusses the potential for a measurement at 
the LHC both with and without the momentum reconstruction technique. Possible 
experimental factors that need to be considered and our final results are explained 
in Section ^. Finally, the Appendices list the Lagrangian, resulting couplings and 
give explicit matrix elements including the full spin correlations that are required in 
our study. 

2. Formalism 

2.1 The process studied and the amplitude squared 

One of the dominant SUSY channels at the LHC is associated squark-gluino produc- 
tion. At the tree level the production process Eq. ( |2.1| ) proceeds via the light quark 
exchange in the direct channel and squark/gluino exchange in the t channel, as seen 
in Fig. I, 

gq Ml- (2.1) 
We study the case where the squark subsequently decays via the following chain: 

Ql ^ + g - Xi^^r + q. (2.2) 

The first step in the cascade decay chain is the two-body process qi — > qX2- Here, 
CP- violating couplings of the X2 enter and are dominated by the phase 0i, see Ap- 
pendix |B.1| . In addition, the spin vector of the X2 has to be explicitly included in 
the amplitude since the full spin correlations have to be taken into account. 

We consider spectra where the second step in the cascade decay chain is the 
three-body decay of the neutralino, X2 ~^ Xi^~^^~ (cf- Appendix p.2|) . The neutralino 
decay occurs via Z° exchange in the direct channel and via slepton i^ ji exchanges 
in the t and u channels, as can be seen in Fig. ^ It is sensitive to CP-violating 



- 5 - 




Figure 1: Feynman diagrams for the associated production process of squarks and gluinos 
at the LHC, gq — > gqi. 




supersymmetric couplings and its structure has been studied in detail in |^ 

cf. Appendix 0. The phase 0i affects the masses of Xi and well as its couplings 

and decay rates. 

Using the narrow width approximation and the formalism of [Q, the squared 
amplitude |Tp of the full process can be factorised into the processes of production 
PP Ml and the subsequent decays — * qX2^ X2 taking into account 

the full spin information of the decaying X2- the purpose of analysing neutralino 
decays we do not need to include decays of the gluino produced with a squark but 
they will be needed later for momentum reconstruction, see Sec. ^ We apply the 
narrow-width approximation for the masses of the intermediate particles, qi and X2j 
which is appropriate since the widths of the respective particles are in all cases much 
smaller than their masses and the mass differences between them are large enough, 
see Appendix 0. The squared amplitude can then be expressed in the form 

\T\' = 4|Afe)|>(x°)pP(m){P(x^)/^(x^) +X^S^p(x^)S"z,(x^)}, (2.3) 

a=l 

where a = 1,2,3 refers to the polarisation states of the neutralino X25 which are 
described by the polarisation vectors ^"(x^)) given in Appendix p.lj . In addition. 
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A(gi) and A{x2) ^^^^ the propagators of the intermediate particles which lead to 
the factors Eq^/niq^Tqj^ and -E-^o/m^oF^o in the narrow-width approximation, 

P{Ml), P{X2) ^{X2) terms in the production and decay that are 

independent of the polarisations of the decaying neutralino, whereas 

Sp(x2) and ^^(xg) are the terms containing the correlations between produc- 
tion and decay spin of the X2- 



According to our choice of the polarisation vectors s°'{x2), see Eq. ( p.4| ) in Ap- 



pendix |RT], Sp/P(x2) is the longitudinal polarisation, Sp/P(x2) is the transverse 
polarisation in the production plane, and Sp/P(x2) is the polarisation perpendicular 
to the reference plane of the neutralino X2- 

2.2 Structure of the T-odd asymmetry 

The tools used in this paper to study CP- violating effects are T-odd observables that 
are based on triple products of the momenta of the involved particles. Here we study 
the visible decay products of the qi in the form of the following triple product that 
can be evaluated in different reference frames: 

^ = Pg- {Pe+ xpe-)- (2.4) 
The T-odd asymmetry is then defined as, 

Nt+ - Nr_ J sign{r}|Tpcilips 



At 



Nr++Nr_ J\T\^d\ips 

J sign{T} ( ELi ^pix'2)^D ix'2))d lips 



J (P{x'2)D{x'2))d\ips 



(2.5) 



where A^^^ (^t_) are the numbers of events for which T is positive (negative) 
and (ilips denotes Lorentz invariant phase space. The denominator in Eq. ( |2.5| ), 
/|Tp(ilips, is equal to the total cross section, namely a{pp qiQ IXi^^^'o)- 
In the corresponding numerator of Eq. ( p.5|) , the triple-product correlations only 
enter via the spin-dependent terms. If the spin of the particles is neglected in the 
calculation the asymmetry will vanish. 

The asymmetry Ar can be visualised as the difference between the number of 
events where the observed Pq (from the decaying qi) lies above (Ar^) or below (A^r_) 
the plane spanned by {p£+,Pe-), normalised by the total number of these events. If 
no complex phases are present, the Pg will on average line up with the plane and no 
asymmetry will be seen. The asymmetry is not a Lorentz invariant quantity but it 
is maximal in the rest frame of the X2- ^t the LHC the momentum of the X2 ^iii 
in general be undetermined. The asymmetry can still be formed in the laboratory 



-7- 



frame but it will be significantly diluted compared to the X2 ^^^^ frame due to the 
initial particle boost that is parametrised by the parton density functions (PDFs). 
The dilution is due to the fact that the relative orientation of the plane spanned by 
Pi+jPi- and the quark momentum can change when boosting from the neutralino 
rest frame to the laboratory frame. Hence, the quark that was initially above the 
plane can instead be observed to come from below the plane in the laboratory frame. 
A possible solution to bypass the effects of the dilution would be to reconstruct the 
momentum of the Xi ^i^nd use it to find the momentum of the X2- ^^^^ discuss 
this method in Sec. ||. 

In order to identify the T-odd contributions, we have to identify those terms in 
|Tp, Eq. ( p.3|) , which contain a triple product of the form shown in Eq. ( |2.4|) . Triple 
products follow from expressions it^ypfja^h^&d'' ^ where a, 6, c, d are 4-momenta or 
spins of the particles involved, which are non-zero only when the momenta are linearly 
independent. The expressions ie^^jpaa^Wc^d'^ are imaginary and when multiplied by 
the imaginary parts of the respective couplings they yield terms that contribute to 
the numerator of Ar, Eq. ( |2.5| ). In our process, T-odd terms with e-tensors are only 
contained in the spin-dependent contributions in the neutralino decay, 

Examining the covariant product, e^^ptja'^b'^cPd" , we can understand why the 
maximal asymmetry is seen in the rest frame of the X2- expand the product 

in terms of the explicit momenta, Eq. (|1.5|) , we see that in general we have triple 
products formed from all different combinations of the momenta of the four different 
particles in the covariant product: 

^i^upaPaPhPcPd = Ea Pb ■ {P^ ^ Pd) + Pd ■ {p^ X pt) (2.6) 

-Eb p^ ■ {pd y< pi) - Ed p^ ■ {pb y- %)■ 

As we evaluate only one triple product we miss the other combinations and the 
asymmetry is not maximal. However, if we are in the rest frame of the X2? 
momentum components of the four vector vanish and we are now only left with the 
single triple product that we are interested in: 



^^lupaPaPbPcPd — ^maPb-ip^xpd)- (2.7) 



Let us finally comment on the origin of the asymmetry, Eq. ( |^.5| ). Neutralinos 
in the process we study are polarised due to the nature of the electroweak coupling 



between the fermion, sfermion and neutralino, see also |2^. Therefore the momentum 
of the outgoing fermion is aligned with the spin of the neutralino X2- "^^^ actual 
asymmetry arises solely in the neutralino X2 decay as a result of the correlation 
between the neutralino spin and the momenta of outgoing leptons, consequently it 
requires us to take into account all spin correlations in the decay chain. 
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Figure 3: The parton level asymmetry At^ Eq. ( |2.5D , for the single decay chain given 
in Eq. (2.2) as a function of the squark momentum, \pq\, in the laboratory frame. The 
scenario is given in Table |l| with the phase (pi = 37r/2. 

3. Momentum reconstruction 



3.1 Dilution effects 

The triple product that is constructed from momenta in the laboratory frame suffers 
from dilution factors at the LHC due to the frame being boosted with respect to the 
rest frame of the ^ more detailed discussion see [Q. It results in a considerable 

reduction in the maximum asymmetry observable when we introduce the PDFs which 
causes an undetermined boost to the system. To illustrate the dilution, Fig. ^ shows 
how the parton level asymmetry Ar-, Eq- ( p.5| ), is diluted in the laboratory frame 
when we produce the q with varying initial momenta. As the X2 is produced in the 
squark decay any boost is transferred. The plot was produced with an analytical 
calculation for the single decay chain given in Eq. (|2.2|) . The scenario displayed in 
Tab. m was used but with the phase set to (pi = 37r/2. 

We see that the asymmetry is maximal. At ~ 14%, when the q is at rest but 
drops to, Ar ~ 2.5%, when \pg\ ~ 1 TeV. The magnitude of the dilution clearly 
depends on the chosen scenario and in particular on the masses of the particles 
involved in the process. If it were possible to reconstruct the momentum of the X25 
we could however perform a Lorentz transformation of all the momenta and bypass 
the dilution factor, potentially recovering the full asymmetry. 

3.2 Reconstruction procedure 

In the mSUGRA scenario that we have chosen to study in this paper (see Section |4.1j ) 
the full reconstruction of the event is made possible by considering the decay chains of 
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Figure 4: The process studied for mo- 
mentum reconstruction. 



both the particles produced in the hard coUision. We include all the particles coming 
from both the ql and the g because there are not enough kinematic constraints to 
perform reconstruction if only the q chain is considered. Therefore we exclusively 
consider the production of ql and g and their subsequent decays, see Fig. ^. The 
(jL decay chain will be the same as considered in Eq. (|2.2|) and the g will decay as 
follows: 

g -*i+t ^ Xt + b + t ^ Xi^+Ui + b + t. (3.1) 

In many scenarios the production of qi along with g is the dominant source of the 
first and second generation squarks at the LHC and in the considered scenarios the 
branching ratios for the decay chain are favourable (cf. Sec 



Assuming that all the masses in the decay chains are known, Sec. [4.1| , the kine- 
matics can be fully reconstructed using the set of invariant conditions and the mea- 
sured missing transverse momentum. For our procedure we follow the methods for 



solving the kinematic equations very closely to those presented in [^5|, The nov- 
elty here is the inclusion of three-body decays of sparticles and allowing for additional 
missing momentum due to neutrinos in the final state. The difference lies also in 
using the mass constraints. For our purpose we assume that masses of sparticles are 
known and aim at the reconstruction of the momenta, whereas the previous studies 
used the above conditions to reconstruct masses. 

Rather than fully reconstructing the kinematics of both decay chains, an alter- 
native idea may be to estimate the momentum of the X2 then boost into this 



approximate frame. A formula that estimates the momentum was presented in 
and is shown below, 

+ (3.2) 
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This approach does not work in our study however as the approximation only becomes 
vahd when p^o — > in the rest frame of the At this kinematical endpoint the two 
leptons are back to back which causes the plane spanned by pi+ and pi- to become 
badly defined. The triple product is therefore small and inaccurate leading to an 
asymmetry that is close to zero. Therefore the approach is not valid in our case. 

In our process the following invariant equations can be formed. 

• q decay chain: 

m|o = (P,oJ^ (3.3) 

^| = (^x?„ + ^.+ +^.„-)'' (3-4) 

m?=(P^o + P,f. (3.5) 

• g decay chain: 

m| = (P^+ + P,)^ (3.7) 

m? = (P^ + PO^ (3.8) 

where P denote 4-momenta of respective particles, and where necessary we 
label particles coming from squark and gluino decays with subscripts a and b, 
respectively. 

• We also have the missing transverse momentum constraint: 

fliss = flo^ + fli, + fl, ■ (3.9) 

• We combine the momenta of Xib ^^^^d coming from the gluino together as it 
is impossible to resolve these two particles: 

P^ME = P^o^ + P., . (3.10) 

• An additional condition on the solutions is that the invariant Pp^g has to be 
greater than the mass of the LSP, m -o : 

P^ME > m\o . (3.11) 

We apply this condition to each solution and discard as unphysical any that 
does not meet it. 
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After expressing the momenta of intermediate particles in terms of the final state 
particles 



Ala 



P,+ + P, 



la ' 



PgME + Pe+ 5 

Pb + P.0 + R+ + P, 



(3.12) 
(3.13) 
(3.14) 



Pi - ^0 I - xi't 

we now have a total of eight equations, Eq. ( |3.3| ) - ( |3.9|) , and eight unknowns: 

EgME,P9ME{x),PgME{y),P9AIE{z)) ■ (3.16) 



In principle the system can be solved to find P^o and PgME- Equations Eq. 
and Eq. (|3.6| ) are quadratic in P^o and PgME respectively, so we consider these last. 
Using on-shell conditions, quadratic terms in the remaining equations can be removed 
giving 6 linear equations, therefore we can simply use a matrix to give us solutions 
in terms of the energies E^o and EgME- We first expand "p^^o and J^gME terms of 
other momenta contained in the respective decay chains: 



' Ala 

l^gME 



(3.17) 
(3.18) 



We can now form the system of 6 linear equations for unknowns A-F with E^o and 
EgME as free parameters 



M 



(A\ 




B 




C 




D 




E 









( h 



X2 



m 



o) + Pp+ ■ Pf- + E^O {Ep+ + En- 

Pa + P,- ■ Pa + E^o Ea 



i(m2 + - m?) + ■ P?, + E^Afs^fc 
|(m? + ml - m]) + P^+ ■ P* + Pj, ■ P* + E^ME^t 



\ 



Pmiss (y) 



(3.19) 



where 



M 



I (tit + ti-)^it i^it + te-)ti- {tit + tt-)t. 












Pi+i^) 



PqPq 





















Pf~ [x) 



Pg{x) 

Pqiy) 



PbPi+ PbPb PbPt 

h 



Pb[X) pt[x) 

Pi+{y) Pb{y) Pt{y) J 

(3.20) 
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The matrix Ai can then be diagonahsed to give solutions for each of the momentum 
components of "jf^o and J^gME terms of E^q and EgME- These solutions are 
substituted into the two quadratic equations, Eq. ( p.3|) and Eq. ( p.6|) , to produce 
two equations of the form: 



Oii-E^o^ + auE^o^EgME + o.22-£^gAf£; + ^^^x'ia ^"^^qME + a = F4 = , (3.21) 
&ii-E^|o^ + buE^o^EgME + b22E'^ME + ^i-^x?a + hEgME + b = Eb = , (3.22) 



where the coefficients aij, ai, a and bij, bi, b are functions only of masses and measured 
momenta. We use, 

Ea-^xEb = (3.23) 

Oil 

to produce the linear equation for E^o , 

anb — abii — a2biiEgME + 0'iib2EgME — CL22bnE?]^^^ + anb22E?j^,j^ , , 

EyO = . (3.24) 

-fliiOi + ctiOii + aubuEgME - anbuEgME 



This result can then be substituted into Eq. ( |3.21|) to obtain a quartic equation of 
the following form: 

Q^^E^ME + QsE^ME + Q'^E'Ime + QlE-gME + Qo = , (3.25) 
where the various Q's are functions of the a's and 6's in Eqs. (|3.21| ) and (|3.22|) . 
3.3 Discussion of graphical solutions 



Analysing the roots of the quartic equations, Eq. (|3.25| ), we select the solutions 



that are real and discard the solutions that contain imaginary parts. The selected 
roots are substituted into Eq. (|3.24|) to find the corresponding solutions for E^o^ . 



Using the values of EgME and E^a together with the diagonahsed matrix we can 
now calculate A, B,C, D, E, F, see Eqs. (|3.17| ),( pn8D and hence the components of 
tgME and t^o^ 



In general, taking into account multiple roots, Eq. ( |3.25|) has 4 solutions. Thus 



we would have 4 real, 2 real and 1 complex pair or 2 complex pairs of roots. Only the 
real roots can yield physical solutions, therefore for each event we normally expect 
real solutions. As Eq. ( p.25| ) is derived from Eqs. ( p.21| ) and ( p.22| ) they share the 



same set of solutions. Both Equations ( |3.21| ) and ( ^.22 ) are polynomials of degree 2 



in E^o and EgME so they correspond to degree 2 curves in the E^o and EgME plane 
(ellipses, hyperbolas or parabolas). The intersection points solve simultaneously both 
equations and at the same time the quartic Equation ( p. 251) . A graphical solution 



to this set of equations is shown in Fig. |^ for one sample event. In this example we 
have two real solutions, of which only one corresponds to the actual momenta in the 
event. 
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Figure 5: Ellipses in the E^o and EgME plane corresponding to Equations ( 3.21] ) (solid 
line) and ( |3.22| ) (dashed line) for one sample event. Out of the two real solutions the lower 
right one gives the correct momenta of the original event. 



In the realistic physical case one has to include uncertainties on measured mo- 
menta and masses. This will of course affect above equations and solutions. Smearing 
of momenta will typically result in shifting of the curves shown in Fig. |^. Therefore 
it is possible that one can get two additional solutions or no real solutions at all. The 



consequences of experimental effects on our analysis will be discussed in Sec. 5.2 



3.4 Practical approach 

In general more than one solution to the above system will remain depending upon 
the resolvent. In some cases Eq. reduces the number of solutions further. For 

the application of this study we therefore calculate the triple product, Eq. (|2.4| ), for 
each individual solution. If all solutions produce the same sign triple product we keep 
the event and use it in our asymmetry. If any of the solutions disagrees on the sign 
of the triple product we will discard the whole event and it will not contribute. This 
method has the disadvantage that we lose events and therefore statistical significance. 
However, we note that it does not introduce any bias, i.e. it can be checked that the 
asymmetry for accepted events remains the same. An alternative method we propose 
is to weight each solution with matrix element and phase space factors to select the 
solutions with the highest probability of being correct. A possible advantage of this 



method is that no events will be lost and this idea is followed up in |^ 



An additional issue when trying to complete the momentum reconstruction pro- 
cedure are combinatorial problems when trying to assign the measured momenta to 
the correct particle in the given event. If we take the leptons as an example, we 



- 14 - 



know that opposite sign, same flavour leptons must come from tlie X2- However, it 
is possible tliat a same flavour lepton could be produced from the xt in the opposite 
decay chain and be confused with those coming from the X2- order to assign the 
leptons correctly to the decay chains one can use the conditions for invariant masses. 
In the squark decay chain a useful observable is the invariant mass of two leptons 
and in the gluino decay chain one can use the invariant mass of the 6-jet from the 
stop decay and the lepton from chargino decay. The end points of these distributions 
are given by 

^e+e- < ^x°o - ' (3-26) 



< 

b 



{ml - ml+ - ml){mlo - m^o) 

(3.27) 



Xl X2 



It turns out that only a small fraction of the events, around 5% in our scenario, sat- 
isfies both conditions simultaneously for the two possible assignments of the leptons 
in question. Moreover, if we run the momentum reconstruction algorithm, we find 
that physically acceptable solutions are only found in roughly 10% of events where 
a wrong assignment has occurred. These conditions act as a strong discriminant, 
therefore we do not expect lepton combinatorics to be a relevant issue for this study. 

Another possible approach to this problem would be the subtraction of the 
opposite-sign opposite-flavour (OSOF) lepton pairs. In that case the quantities Nj-^ 
and Nr_ in Eq. (|2.5| ) would be defined as 

iVr+ = [iVe+e- + N^+^- - N,+^- - N,-^+]r^ , (3.28) 
Nr. = [Ne+e- + N^+f.- - N,+^^ - N,^^+]r_ . (3.29) 



In Eqs. (|3.28|) and (|3.29| ) one adds the number of events with a positive (negative) 



triple product using e'^e~ pairs (and the jet) to the number of events with a positive 
(negative) triple product using pairs and subtract the number of events with 

positive triple product using ^~ and e~/i'^ pairs. For the combinatorial background 
one gets equal rates for same flavour and opposite flavour lepton pairs since the 
leptons are uncorrelated. Hence, this "flavour subtraction" procedure removes (up 
to statistical fluctuations) all combinatorial background resulting from the lepton 
pairing problem and also all background coming from X2 ~^ Xi^^t^- 



4. Numerical results 

4.1 Chosen scenario: spectrum and decay modes 

In order to study the experimental prospects of measuring CP-violating effects at 
the LHC we use an MSSM scenario derived from mSUGRA parameters defined at 
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Parameter 


Value 


Particle 


Mass [GeV] 


Particle 


Mass [GeV] 


mo 


150 GeV 


9 


496.5 


X? 


78.1 




200 GeV 


di 


484.1 




148.4 


Ao 


-650 GeV 


d-R 


466.4 


xf 


148.2 


tan (3 


10 


ul 


477.9 


xt 


436.0 


sign /i 


+ 


ur 


465.9 


ei 


207.5 


\Mi\ 


80.5 GeV 




397.2 


eR 


173.1 


M2 


153.3 GeV 


k 


462.6 


i>e 


192.0 


Ms 


484.6 GeV 


ii 


171.0 


ri 


149.4 




419.0 GeV 




498.0 


T2 


212.5 
187.2 



Table 1: mSUGRA input parameters at the GUT scale, MSSM parameters and particle 



masses in GeV from SPheno 2 . 2 . 3 |4S] for = and with mt = 171 GeV. 



the GUT scale, as shown in Tab. |l|. This scenario has been already used to analyse 



the properties of the neutralino sector in [31]. The values of the parameters at the 



electroweak scale have been derived using the RGE code SPheno [^. Masses of 
the coloured particles are of order 500 GeV, apart from the hght stop ti, which 
has a mass of 171.0 GeV. The lightest supersymmetric particle (LSP) is the lightest 
neutralino with a mass of 78.1 GeV. The second neutralino and the light chargino 
have masses around 150 GeV, whereas the sleptons are around 200 GeV. Therefore 
we note that the two-body decay channels of X2 ^'^d xf are closed and this gives 
a good opportunity to study CP-violation effects in their three-body decay modes. 
Details of the mass spectrum can be found in Tab. |l|. The values of the gaugino mass 
parameters reproduce the given spectrum in the case when all the CP phases are set 
to 0. In order to generate CP-violating effects in the following we will attribute a 
non-zero phase 0i to the bino mass parameter, 

Mi = \Mi\e''^\ O<0i<27r, (4.1) 

while keeping the absolute value \Mi\ fixed as given in Tab. 

Although we chose a specific scenario our method is applicable in a wider range of 
parameter points. For the decay chain that we concentrate on, we require the three- 
body decay of the X2 which places the following constraints on the SUSY masses, 

m^o < mi , (4.2) 
m^o — mj^o < mz- (4.3) 
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Equation (|4.2|) ensures that the decay, X2 ~^ cannot occur whereas Eq. (|4.3| ) 
forbids the decay X2 Zxi- 

To complete momentum reconstruction from the gluino side, other decay modes 
can contribute with regard to the scenario presented here but their impact on the 
final result would be scenario dependent. For example if the decay t became 
kinematically open the momentum reconstruction would actually become easier as 
there would be no neutrino in the final state and the system would become over- 
constrained. However, there would now be one additional top and one additional 
lepton in the final state that may cause new combinatorical difficulties. 

As already mentioned, the highest production rate is typically obtained for 
coloured final states containing squarks and gluinos at the LHC In our scenario, 
where their masses are not very heavy, the total cross section for production of 
strongly interacting supersymmetric particles reaches 140 pb, see Tab. H. For our 
purpose we will be interested in the inclusive production of left squarks and the 
associated production of left squark and gluino^. As was mentioned in Sec. |I], anti- 
squarks give a CP asymmetry with exactly the opposite sign to squarks. We note 
however that the inclusive cross section for left squark production is almost a factor 
of 4 larger than the cross section for left anti-squark production. This is a direct 
consequence of the fact that we have two protons in the initial state for which the 
abundance of quarks is significantly higher than of ant i- quarks. Since we cannot 
distinguish experimentally squarks and anti-squarks of the first two generations this 
fraction of anti-squarks will cause some dilution in the asymmetry that will be how- 
ever compensated for by a very high production rate. A similar situation occurs for 
associated squark gluino production for which the ratio of squarks to anti-squarks is 
18.2 pb to 3.1 pb, see Tab. 0. Together with left squarks we also have the produc- 
tion of right squarks at a comparable rate. However, since the latter decay almost 
exclusively to the lightest neutralino in our scenario, as shown in Tab. ^ they do not 
give rise to the CP-odd asymmetry. 

Following the production process we include subsequent decays of squarks and 
gluinos. The dominant decay mode of the gluino is to the light stop and the top with a 
branching ratio of BR = 53.8%. The light stop then decays almost exclusively to Xi^- 
Left squarks decay mainly to the light chargino and the second lightest neutralino 
with branching fractions of 65% and 33% respectively. Finally we consider the decays 
of chargino Xi ^'^^ neutralino X2- Leptonic decays constitute in total BR = 61% 
for the chargino and BR = 68% for the neutralino decay modes. For decays into 
light leptons we have BR = 24.3% and BR = 9% for chargino and neutralino, 
respectively. The most interesting are the neutralino decays to electrons and muons 

^Triple product asymmetry in right squarks decay chain has the opposite sign as compared to the 
left squarks (and the same sign as that of left anti-squarks), however in our case where the coupling 
of the right squark to the second lightest neutralino is suppressed due to small bino component of 
X2, this contribution remains negligible. 
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Produced Particles 


Cross Section (pb) 


At least one coloured SUSY particle. 


148 


At least one g. 


58.8 


qlQl + qiql + Mr + (iLq*R + qig 


30.0 


qlqL + qlq*L + qlqn + qlq*R + qlg 


8.3 


qig 


18.2 


ql9 


3.1 




66.3 



Table 2: Leading order cross sections at -y/s = 14TeV for direct production of various 



particles from Herwig++ 2.3.2 |4^, |4l[ using MRST 2004L0 PDF set |50[. q stands for 
squarks of the first and second generation. 



Decay 



g^tit + c.c. 
g — > hib + c.c. 
g + c.c. 
g -^b2b + c.c. 
g qLq + c.c. 



qR 
qR 



Xiq 
xlq 



BR 



53.8% 

26.6% 

11.8% 

3.8% 

3.3% 



98% 
1% 



Decay 



qL 

qi 
qi 



Xi 

xt 

xt- 

xt 



xfq 
X2q 
x\q 



' x\quqd 



BR 



65% 
33% 
1.5% 



37.2% 
38.5% 
12.2% 
12.1% 



Decay 



ti ^ xt^ 
ii ^ 



X'2 



xl 
xl 



xl 



xlc 



xl-r-T^ 

x\qq 



BR 



98.1% 
1.6% 



59.3% 
23.6% 
8.1% 
4.5% 
4.5% 



Table 3: Branching ratios for the scenario defined in Tab. | from SPheno 2.2.3 Q for 
^1 =0. 



that are used to construct the CP-sensitive triple product, Eq. (p.4|). The summary 
of the relevant branching ratios can be found in Tab. |^. It may be noted that in 
the scenario presented, the branching ratio X2 ~^ Xi'^^^''"" is unusually high (59.3%). 
This is due to the masses of the X2 ^ being very close and hence the kinematic 
factors are favourable. The particular branching ratio is simply a coincidence in the 
scenario chosen however and is not a required feature for our study. 

When we vary the phase of the Mi parameter, the masses and couplings in 
the neutralino sector are affected. First we note that the changes in the neutralino 
masses are negligible and smaller than the possible experimental accuracy. It turns 
out however, that the phase 0i has large impact on neutralino couplings and therefore 
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the pattern of its decay modes. The most significant change is for the hght chargino 
decays to the LSP and a fermion pair. With increasing phase the branching ratio 
for xt ^ Xi^tT^ rises and eventually reaches 70% for (pi = vr. At the same time 
the branching ratios for decays to light leptons remain roughly at the level of 10%. 
As we will explain later, a chargino decay to tau, followed by a leptonic tau decay 
can be used for momentum reconstruction in the same way as a chargino decay to 
an electron and a muon. On the other hand, the decay X2 ~^ Xi'^~^'^~ followed by 
a leptonic tau decay will be a background, since it gives incorrectly reconstructed 
momenta. Finally, we note that in the case of neutralino X2 decays to light leptons, 
the respective branching ratios increase up to 5.5% for 0i = 7r/2. 

Since our procedure requires the sparticle masses to be known, we wish to com- 
ment on the possibility of the mass determination in our scenario. One of the stan- 
dard approaches at the LHC is to study kinematic edges and endpoints of the invari- 
ant mass distributions, see e.g. In our case the masses of g, ql, h, X12 xt 
are required. The possible invariant mass observables would include: 

• quark and leptons in the decay chain of followed by the squark decay to 
chargino or neutralino 



lepton pair in the neutralino X2 decay, cf. Eq. ( p.26[ ) 



top and bottom quarks in the gluino decay chain Eq. (|3.1| ), 
bottom quark and lepton in the stop decay chain, cf. Eq. (|3.27|) , 



• quark pairs in the decay of gluino to the 1st and 2nd generation squarks |]52| . 

Fitting the above invariant masses to experimentally measured edges and endpoints 
should provide enough number of constraints to obtain required information. This 
method can be perhaps supplemented by a mass reconstruction method, see for 
example p^, |38|, B^, EOl HQ, EM. Therefore, in the following we assume that 



the masses are known. The impact of the mass uncertainties will be discussed in 
Section 



4.2 CP asymmetry at the parton level 

We start by discussing the dependence of the asymmetry on 0i, Eq. (p.5|), at the 
partonic level studying only the decay chain presented in Eq. ( p^.2| ). In order to see 
a maximal effect, we place the at rest and calculate the triple product and the 



asymmetry in its distribution. Sec. ^ 



We see from Fig. ^ that the asymmetry in this scenario is roughly 14% and this 
occurs when the phase 0i is just above it/ 2 and just below 3tt/2. The asymmetry 
is produced by a complex interplay between different couplings in the X2 decay, cf. 
Eqs. ( p.24|) - (p.26| ), and can vary significantly between scenarios. These couplings can 
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Figure 6: The parton level asymmetry At, Eq. ( |2.5| ), for the single decay chain given in 
Eq. ( |2.2| ) in the rest frame of qi as a function of The scenario is given in Tab |l|. 



all have different behaviour with respect to 0i and in other scenarios the maximum 
asymmetry can be seen far from 7r/2. 

It should be noted from the previous plot that the asymmetry is obviously a 
CP-odd quantity that in addition to a measurement of the phase, also determines 
whether it lies above or below vr, as seen in Fig. ^. In comparison, using CP-even 
quantities, for example the mass, it is not possible to determine if the phase is above 
or below vr. It must also be noted that these quantities have a weak dependence upon 
the phase and will be challenging to study. Moreover CP-even quantities alone do 
not give unambiguous proof of CP violation in the model, that can only be provided 
by CP-odd observables. 



4.3 Influence of parton distribution functions on CP asymmetries 

Experimentally the situation becomes significantly more complicated since in general, 
particles are not produced at rest but can be heavily boosted in the laboratory frame. 
Our observables can be significantly reduced in size by this effect as triple product 
correlations induced by spin effects are largest in the rest frame of the decaying 
particle. Essentially, a boosted frame can make the momentum vector of the quark 
appear to come from the opposite side of the plane formed by and As explained 



in Sec. p.2| this causes a severe dilution in the asymmetry that is measured at the 
LHC. 

There are two other dilution factors that have to be taken into account and 
give a further reduction in the observed asymmetry. The first one we consider is 
the contamination from anti-squarks g2 that will be produced at the LHC along 
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with squarks q^. There is no way of identifying the charge of the decaying particle 
and consequently we need to include the effect of the g2 in our analysis. Anti- 
squarks produce an asymmetry of the opposite sign to qi so if for example we had 
equal numbers of each of them, no overall asymmetry could be seen. However, the 
production cross section for g2 is substantially lower than that of qi due to the valence 
quarks in the proton, see Tab. |^ and at the LHC we would expect only roughly 20% 
of the sample to be g2. 

The other background contribution is that of neutralinos decaying to a tau pair, 
X2 ~^ Xi^^'T^ followed by leptonic tau decays to the pair of opposite-sign same- 
flavour leptons. Since the branching ratio for the above decay is large, even after 
inclusion of leptonic tau decays there is a significant number of events faking our 
signal, i.e. 

BR{xl xlr+T- x^^i^rueUiUrUr) ^ 2 X 0.6 X 0.175 x 0.175 = 3.7%. (4.4) 
compared with 

BR{xl ^ XiCj:,,) ^ 2 X 0.045 = 9%. (4.5) 

As the asymmetry calculated for such leptons is diluted, this background introduces 
a further dilution factor. However, this background can be easily removed using the 
flavour subtraction technique described above, Eqs. (|3.28|) and (|3.29|) . In addition 
we will see later that much of this background is removed after the application of 
simple selection cuts on the lepton energy and the invariant mass. 

We use the MRST 2004L0 |j5^ PDF set in our analysis of the asymmetry and 
plot the integrated asymmetry At as a function of 0i in Fig. 0. We see that the 
inclusion of the PDFs and the sample reduce the asymmetry by about a factor 
of 8 in this scenario as compared to the result of Sec. |4.2| . The maximum asymmetry 
is now \At\ = 1.7%. It must be noted that the dilution factor does depend on 
the scenario studied and changes in particular with the mass of the sparticle that is 
initially produced. 

For the calculation of the asymmetry we included the production channels shown 
in rows 3 and 4 of Table |^ but we only take decays of individual qi and g2 following 
the decay chain in Eq. ( |2.2|) . At this point correct identification was assumed for the 
final state particles and no hadronisation or detector effects were included but these 
assumptions will be relaxed in Sec. ^ The only backgrounds in the study are those 
discussed above. 

Using the total production cross sections^ for (Jl and g2, and respective branching 
ratios from Table |^ we can now estimate the integrated luminosity required to observe 

^Note that the total rate of squark production is actually larger than 30.0 pb quoted in Table 
as some of the subprocesses give a pair of squarks e.g.{qLqL) and both can contribute to our analysis. 
The total number of squarks in the sample at the given luminosity C and ^/s = 14 TeV is therefore 
33.2 X C. In case of anti-squarks the number is 8.5 x C. 
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Figure 7 : The parton level asymmetry At, Eq. (|2.5| ), in the laboratory frame with PDFs 
included in the analysis using the scenario shown in Tab. |l[ The coloured lines show the 
size of the asymmetry needed for a 3o" observation at the given luminosity, £=(20 fb"\ 50 
fb"^, 100 fb~^), assuming squarks were produced via the channels shown in rows 3 and 4 
of Tab. |2|. All produced q and q* that follow the decay chain given in Eq. ( |2.2| ) are taken 
into account. 




the asymmetry at the LHC with a certain significance. We assume that {Nq-__), 
the numbers of events where T is positive (negative) as in Eq. ( p.5|) , are binomially 
distributed, giving the following statistical error ||55| : 

^(^^)stat ^ 2^/eil-e)/N, (4.6) 

where e = Nr+/{Nr+ + Nr_) = \{l + At), and = Nr+ + A^r_ is the total number 
of events. This can be rearranged to give the required number of events for a de- 
sired significance. The horizontal lines in Fig. |^ show an estimate of the amount of 
luminosity required for a Scr observation of a non-zero asymmetry. In other words, 
an asymmetry can be seen at the Sex level where the asymmetry curve in Fig. |^ lies 
outside the luminosity band. 

4.4 Impact of momentum reconstruction on the observable CP asymmetry 

In order to increase the statistical significance of our CP asymmetry we investigate 
the possibility of reconstructing the momenta of the invisible particles in our process. 
In principle a perfect reconstruction would return the magnitude of the asymmetry 
to that where the qi is at rest but in reality there are additional complications with 
this procedure, see Sec. The reconstruction is performed at the partonic level with 
PDF's included in the production process. 
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Figure 8: The parton level asymmetry At, Eq. ( |2.5| ), in the reconstructed qi rest frame. 
The coloured lines show the size of the asymmetry needed for a 3a observation at the given 
luminosity, £=(20 fb^^, 50 fb^^, 100 fb^^) for the production channels qig, q^g Tab. |2|. 



The decay chains included are shown in Eq. ( |2.2| ) and Eq (3.1). As explained in the text 
correct jet and lepton assignment is assumed. 



Again our sample of events will contain q1 which have an asymmetry of the op- 
posite sign to that of qi as has already been discussed in Sec. ^j3[ As we are looking 
exclusively at qig and qig production however when applying the momentum recon- 
struction, we actually have a smaller number of that dilute the asymmetry (15%, 
see Table H). After including this dilution factor we see our maximum asymmetry 
reduced from Ar ~ 14% to At ~ 11%, see Figs. H and H, respectively. 

To calculate the luminosity we require to see a statistically significant effect at 
the LHC we include the production cross sections for channels qig, q^g (Tab. 0) and 
the branching ratios from both the qi and g decay chains Eq. ( p. 21 ) and Eq. ( p.lj ). No 
hadronisation or detector effects were included in this section and correct identifica- 
tion was assumed for the quarks. For the leptons the correct assignment is made as is 
explained in Sec. ^]4[ For the jets, a detailed experimental study would be required 
to look at this question as we need to include hard radiation, reconstruction and 
6-tagging efficiencies. However, a cursory examination suggests that the individual 
jets may be resolvable. For example, the jet coming from the qi decay is the hardest 
one for a high proportion of events. We would also require three other jets to have 
the invariant mass m^. Finally we would then require at least one 6-jet to be tagged. 
We would like to state here that the only backgrounds included in the study are 
those of q]^g production and the decay chain Eq. (|4^ ). 

After inclusion of all branching ratios the production rate for our process drops 
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down to 200 fb (after hadronic top decay), which results in approximately 20 000 
events at the integrated luminosity L = 100 fb~^. One extra point to note is that 
we are able to use events where the xf the opposite decay chain produces a 
r^. Kinematically these events are similar to our normal signal events but we have 
one extra neutrino u^- that will be invisible and simply contribute to our missing 
transverse momenta. An additional factor that reduces the number of reconstructed 
and accepted events are the multiple solutions, as discussed in Sec. ^. We only use 
events when all solutions produce the same sign triple product and this allows us to 
use approximately ~ 60% of the events. 

The horizontal lines in Fig. |^ again show an estimate of the amount of luminosity 
required for a Sex observation of a non-zero asymmetry. It can be seen that the 
luminosity lines are not flat but actually give a narrower band as (pi vr. This effect 
is caused by the branching ratios, — > qX2^ X2 ~^ Xi^~^^~ xt ~^ Xi^d~ altering 
with the change in phase and producing more events. We see that after 20 fb~^ of 
well understood LHC data it may be possible to start seeing a statistically signiflcant 
effect if large phases are present. It must be noted however that for this method to 
be successful we will require mass measurements of the individual particles involved 
in the decay chains. This will obviously require signiflcant running time and may 
even need the help of a precision linear collider. We will address this issue in more 
detail in the next section. 



5. Inclusion of experimental factors 

To more realistically estimate if the study will be possible at the LHC some additional 
experimental factors need to be included in the analysis. For this purpose we simulate 
events using the Monte Carlo event generator Herwig++ 2.3.2, that has all the spin 
correlations included as required in our analysis. These events will be used in the 
following to perform momentum reconstruction and as a cross check for our analytic 
calculations. 

We start with the inclusion of the selection cuts that have to be used to resolve 
leptons and jets and these are listed below: 

Et{ji) > 100 GeV, (5.1) 

^t(j)> 25 GeV, (5.2) 

^t(4,^) > 10 GeV, (5.3) 

Me+e^ > 20 GeV , (5.4) 

\v\ < 2.5. (5.5) 

Here Ej-iji) is the transverse energy of the hardest jet, Ej-ij) is the transverse energy 
of all other jets, ET{ie,^i) is the transverse energy of the leptons, Mi+i- is the invariant 
mass of the opposite-sign same flavour lepton pair, and 77 is the pseudorapidity of all 
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the final state particles in the decay chain. Moreover we require at least two 6-jets 
and that 1 6-jet plus 2 other jets (typically with the lowest pt) should reconstruct 
the top quark. Since we need the top momentum to be reconstructed we only take 
into account tops that decay hadronically. 

One of the consequences of the application of the above cuts, especially Eq. ( |5.3| ) 
and ( |5.4| ) is the significant reduction in the background originating from r's, Eq. ( [4.4| ). 
This is due to the rather low energy of the leptons coming from r decays and the 
even lower invariant mass of the resulting lepton pair, which is peaking around 0. 
Already at this point, approximately 95% of this background is removed. 

Another factor we include is the momenta of the resolved particles being smeared 
due to the intrinsic experimental precision. The accuracy for both jets and electrons 



follows the same function but with different coefficients 56 



' ^ + + 4e , (5.6) 



E y E ' ' "^'^ ' 



where 



. 1 



for jets ttj = 0.6 GeV2, bj = 1.5 GeV and Cj = 0.03; 



1 



• for electrons the accuracy is better, with = 0.12 GeVa, be = 0.2 GeV and 
Ce = 0.01. 

The resolution for muons has a different functional dependence, 

cTp^ r 0.00008(pT/GeV - 100) + 0.03, pr > 100 GeV, 



Pt { 0.03, pT<100GeV. ^^''^'^ 

In addition we also have a finite resolution on the measurement of missing transverse 
energy, 

Et Et ^Er/GeV ' 
which are the errors on the x and y components of the MET vector and Et is the 
scalar sum of all visible transverse energy . 

The momenta smearing will only affect the observables when we perform momen- 
tum reconstruction as the LSP momenta will be reconstructed with limited precision. 
The triple product will not suffer however as the measurement only relies on the di- 
rection of the measured particles and not on the energy measurement. The direction 
can be found far more accurately and this error happens to be negligible for our 
observables. 

Finally, we also investigate the fact that the masses of the particles in the decay 
chains we are interested in will only be known with a certain precision at the LHC (we 
assume 10% error), see Section ^TT]. This error will again only affect the observables 
when we perform momentum reconstruction in order to boost into the rest frame of 
the X2 and can cause the frame to be mis-measured, see Section p.2| for more details. 
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Figure 9: The asymmetry Ar, Eq. (p.5|), in the laboratory frame for the decay chain 
Eq. ( |2.2| ) after the cuts, Eqs. (^?l|)-( [5^ ), have been apphed. The coloured lines show the 
size of the asymmetry needed for a 3cr observation at the given luminosity, £=(50 fb^^, 
100 fb-\ 300 fb-i), assuming squarks were produced via the channels shown in rows 3 
and 4 of Tab. ^. Momentum smearing for both the leptons and quarks was studied and 
found to have a negligible effect. All other relevant experimental details are mentioned in 
the text. 




5.1 Experimental factors without momentum reconstruction 

Out of the experimental factors mentioned above, only the cuts affect the result for 
the triple product correlation measured in the laboratory frame. These cuts reduce 
the number of detectable events by ~ 50% and consequently significantly increase 
the luminosity required to make a statistically significant measurement at the LHC, 
see Fig. For example if large phases are present we may begin to see hints with 
integrated luminosity £ = 50 fb~^. With C = 300 fb~^ we could become sensitive 
to phases in the ranges 0.15 tt ^ 0i ^ 0.9 vr and 1.1 vr ^ 0i ^ 1.85 vr where the 
asymmetry \Ar\ > 0.7%. 

For the calculation of the asymmetry we included the production channels shown 
in rows 3 and 4 of Table ^ taking into account the decays of individual ql and g2 as 
listed in Eq (|2.2|). No explicit hadronisation was included but momentum smearing is 



expected to simulate some of this effect. We also need to include additional hard QCD 
radiation and other detector effects (for example fakes) in a full experimental study. 
Correct identification was assumed for the quark and two leptons. The backgrounds 
in the study are those from and Eq. ( ^4.4|) . 

It must be noted however, that significant pollution due to backgrounds will 
be expected for this signal from both the Standard Model and the MSSM. Further 
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experimental cuts will certainly be required to improve the signal/background ratio 
and without a more detailed study it is hard to predict what effect this will have on 
the ability to complete this measurement at the LHC, especially as the asymmetries 
are rather small. However, a more detailed experimental analysis is beyond the scope 



of the present paper [HS 



5.2 Experimental factors with momentum reconstruction 

When we perform the momentum reconstruction we need to include the experimen- 
tal precision on the momentum of the visible particles. This resolution is ~ 3% 
for leptons and follows Eq. ( ^.6|) for jets. The corresponding effect on momentum 
reconstruction is a reduction in the number of events that have the same sign triple 
product. As stated in Sec. [4.4| we discard any events where we have solutions with 
opposite sign triple products. Discarding these events reduces the percentage we can 
use from ~ 60% without momentum smearing down to ~ 30% when we include it. 

The other difficulty momentum smearing creates is that all the reconstructed 
solutions can now have the wrong sign triple product as we no longer correctly 
reproduce the rest frame of the neutralino X2- Inevitably this effect produces a 
decrease in the observed asymmetry from ~ 11% to ~ 8%. 

We again include the cuts on all visible particles in our decay chain given by 
Eqs. ( ^.1|) -(|575| ) . These cuts significantly reduce the number of visible events and 
remove ~ 80% of the events compared with our initial naive estimates. When we 
combine the cuts with the momentum reconstruction efficiency we are left with ~ 6% 
of the initial events and this clearly increases the luminosity needed to make an 
observation statistically significant. After inclusion of these effects the number of 
events drops from 20000 down to 1200 at the integrated luminosity of £ = 100 fb~^. 
This results in a la absolute uncertainty of order ~ 3% on the asymmetry, according 
to Eq. dJ). 

Another possible experimental aspect we investigate is a 10% uncertainty on the 
masses of the supersymmetric particles used in the momentum reconstruction. We 
found that this has a negligible effect on the momentum reconstruction as long as the 
mass differences between different particles in the decay chain are known better than 
(9(5 GeV). The assumption that the mass differences are known to a higher accuracy 
is reasonable as the main method of measuring masses in SUSY decay chains at the 
LHC will be via kinematic end points that are measured with high precision. 

Figure |TU| shows the asymmetry and the luminosity required at the LHC to 
see a statistically significant result at the 3a level at the LHC once all the above 
factors have been taken into account. The production channels are g^^, q}^g and 
the branching ratios are included from both the and g decay chains. No explicit 
hadronisation effects were included here and possible additional hard jets in the final 
state require further study. Also correct identification was assumed for all the leptons 



and quarks in both decay chains as explained in Sec. 4.4 and Sec. p.4|. Again the 
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Figure 10: The asymmetry At, Eq. ( |2.5| ), in the reconstructed qi frame for qLd, 0*^9 pro- 
duction followed by the decay chains shown in Eq. ( p.2| ) and Eq (|3.1|) with cuts Eqs. (5.1)- 
( |5.5| ) . The coloured lines show the size of the asymmetry needed for a 2>a observation at the 
given luminosity, £=(50 fb~^, 100 fb~^, 300 fb^^). The momenta of the final state particles 
have been smeared according to Eq. ( |5.6D -(^^) to replicate the LHC detector effects. All 
other relevant experimental details are mentioned in the text. 



only backgrounds included in the study are those of q}]) production and of taus that 
decay to visible leptons, Eq. ( [4.4| ). 

We can see that using this method, with integrated luminosity C ~ 100 fb^^ we 
start to become sensitive if large phases happen to be present. After C ~ 300 fb^^ we 
expect to have sensitivity to phases in the ranges 0.2 tt < 0i ^ 0.85 tt and 1.15 vr < 
01 < 1.8 TT, and obviously more luminosity will improve this further. 

Note that a direct comparison between the methods with and without momentum 
reconstruction based on the above plots should not be performed. The backgrounds 
from both SM and MSSM will be more severe when we do not perform momentum 
reconstruction and clearly many new cuts will be required to isolate the signal, which 
is the consequence of being totally inclusive. For the method when we perform 
momentum reconstruction, we have a well defined final state that is difficult to 
be faked by Standard Model processes and is also uncommon for SUSY cascade 
decays. Moreover the multiple cuts on all the particles and the need for the missing 
momentum to be successfully reconstructed mean that many backgrounds will be 
rejected. 

5.3 Ongoing work 

A new approach in the future study is to use matrix element weighting for 
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each of the multiple solutions. This would significantly increase the number of ac- 
cepted events and possibly increase the statistical significance of the result as we 
had observed that some of the "wrong" solutions obey rather unlikely kinematic 
configurations. However, the validity of this method needs further study. 

The other possible improvement is the application of the cuts on anti-squarks. In 
the associated production, Eq. anti-squarks come always from the sea quarks, 

whereas squarks can originate both from the sea quarks and the valence quarks. 
Therefore one can expect some difference in their distributions. This would give a 
handle to cut out one of our most severe backgrounds. Prehminary studies showed 
that indeed this selection criteria can turn out to be useful depending on the actual 
scenario and the collected luminosity. 

It is clear that to fully understand if these measurements are possible at the 
LHC a detector simulation including backgrounds will be required. This study has 
to include the combinatorial issues that may be encountered with both the final state 
leptons and the jets. In addition, the impact of particle reconstruction inefficiencies, 
falsely identified particles ("fakes"), different detector acceptance for different par- 
ticle species and other experimental considerations need to be considered. Such an 
analysis is in the line of future work. 

Further improvements in measuring the CP-violating phase (pi would require the 
inclusion of more observables sensitive to the phase, like lepton opening angle or 
lepton angular distribution ||31[. Nevertheless, in order to determine the phase in a 
precise and rather model-independent way, experimental input from the clean and 
unbiased measurements at a future linear collider |3^, |5^, ^ is required in addition 
to the LHC. 



6. Conclusions 

In this paper we have investigated the problem of discovering CP-violating effects at 
the Large Hadron Collider. Our tool in this analysis was the triple product correlation 
of momenta of the quark and the lepton pair from the X2 decay. This kind of T-odd 
observables can give us a handle to discover the presence of complex phases in the 
model. We took into account the full production and decay process of squark and 
gluino. Since triple products depend crucially on spin correlations, they have been 
included both in the analytical calculation and the event generation, that has been 
performed using Herwig++ 2.3.2. The process of special interest in our case was 
the squark decay into quark and neutralino X2 followed by the three-body leptonic 
neutralino decay. For this decay in our scenario one can expect an asymmetry in 
the triple product distribution of up to 14% when calculated in the rest frame of the 
decaying neutralino. The source of the CP violation in our case was the phase of the 
bino mass parameter Mi. 
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Due to the hadronic experimental environment of the LHC, precise measurements 
will be a challenge both from experimental and theoretical point of view. Therefore 
we have considered two possible methods of measuring CP-violating effects in the 
neutrahno decay. The first method calculates the triple product in the laboratory 
frame. It turns out however that the initial CP-odd asymmetry is diluted both by 
high boosts of produced particles and by the admixture of anti-squarks. This makes 
its observation far more difficult but can give us hints of CP violation after basic 
selection cuts in our example. 

As mentioned before the largest effects of CP violation can be expected in the rest 
frame of the decaying neutralino. We studied the use of momentum reconstruction 
of invisible LSPs and neutrinos to get access to the rest frame of the neutralino 
X2- Using a set of invariant kinematic conditions we showed that it is possible to 
fully reconstruct the production and decay process on an event-by-event basis. The 
important result here is the possibility of reconstruction in the case of three-body 
decays and escaping neutrinos. 

Having fully reconstructed events we are able to boost particle momenta back 
to the rest frame of the neutralino X2 the asymmetry becomes as large as 10%. 
Again, we include experimental event selection criteria as well as momentum smear- 
ing to account for finite detector resolution. We also consider the most severe back- 
ground due to anti-squarks and we showed that the signal contamination from neu- 
tralino decays to leptonically decaying tau pairs does not pose a problem for our 
measurement. 

Applying momentum reconstruction and taking into account the experimental 
effects we have found that the CP-odd asymmetry drops down to 8%. To avoid 
ambiguities and due to experimental constraints many events have to be discarded, 
nevertheless the left-over sample is enough to probe CP-violation effects. One should 
see a 3a effect at £ = 300 fb~^ and we expect sensitivity to the phase in the range 
0.2 TT ^ (pi ^ 0.85 TT. We emphasise that the asymmetry after momentum reconstruc- 
tion is a much cleaner observable from a theoretical point of view, thanks to a well 
defined final state. Therefore, using the above technique there are good prospects 
for the observation or exclusion of CP- violating effects for a wide range of the phase 
01 after a few years of LHC running at the high luminosity. The full assessment of 
LHC's ability to resolve CP violation in the MSSM will definitely require a realistic 
simulation of detector effects, SM and SUSY backgrounds. That issue is beyond the 
scope of the present phenomenological analysis and therefore we leave this question 
for future study. 

The momentum reconstruction procedure turns out to be extremely useful method 
for studying new physics models at the LHC, as we have shown in this paper. Nev- 
ertheless, in order to reveal the underlying model, including CP-violating phases, a 
more precise machine, like an e~^e~ linear collider, will be required. 
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Appendices 

A. Interaction Lagrangian and neutralino mixing 

Neglecting the mixing effects in the slepton sector, the interaction Lagrangian terms 
for X2 decay are 

jOzoi+i- = ^—Z/Y[LePL + RiPR]i , (A.l) 

cos 6w 

Z COS u\Y 

>C«^o = gfi, iPn xl h + 9fR^ IPl xl 4 + h.c. , (A.3) 

where Pl,r = |(1=F75), 9 = e/ sin6'iy and 9w is the weak mixing angle. The 
couplings are given by 

fi, = ^itan9wNki + Nk2), (A.4) 

fj,, = -V2t^n9wN;,, (A.5) 

0^^ = N^,N:,-Nm3K,, (A.6) 

= —O^* (A 7) 

= -i + sin^^w', (A.8) 

i?£ = sin^ Ow , (A.9) 

where N is the neutralino mixing matrix defined by Eq. ( |A.13| ). 
The interaction Lagrangian for decay is 

^^qqxl = aflk^PR xl Ql + h.c, (A. 10) 

with 

flk = ^ {~T,,tanewNk, + SiV^s) , (A.ll) 

where T^g is a weak isospin of a squark and we have neglected small Yukawa couplings 

for the first two generations of squarks and mixing between left and right squark 
states. 

The complex symmetric neutralino mass matrix Ai^ is given in the basis {B, W^, H^, H2) 

by: 

f Ml —mzcpsw rnzspsw \ 

M2 mzcpcw -mzSfsCw 

-mzcpsw mzcpcw -/i 
y mzspsw -mzspcw -/U ^ J 



M 



N 



(A.12) 
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where the abbreviations sw = sinOw, cw = cosOw, sp = sin/?, cp = cos/? have been 
used and tan/5 = V2/V1 is the ratio of the vacuum expectation values of the Higgs 
fields. It is diagonalised by the unitary 4x4 matrix 

diag(m-j^o, m^o, m^o, m^o) = N*A4nN'^ , ("^x? ^ '"^x" ^ ^ ^x'O (A. 13) 

giving 4 non-negative Majorana neutralino mass eigenstates. The masses in Eq. ( |A.13| ) 
can be chosen real and positive by a suitable definition of the unitary matrix A^. For 
convenience, in the following we use the notation: rrii = m^o. 



B. Amplitude squared including full spin correlations 

B.l Neutralino production qi — > 

The analytic expression for the production density matrix can be decomposed as 

\M{qL - XW = Pifjq) + S?>(X?) , (B.l) 
whose spin-independent contribution reads 

P{x1q) = ^-^\fl/{p,P^o), (B.2) 

where Pg and p^o denote the four-momenta of the quark q and the neutralino Xj- 
The spin-dependent contributions is T-even and given by 

S?>(x;) = -yl/l/m,(p,.»(x?)), (B.3) 

where s"(Xj) denotes the spin-basis vector of the neutralino Xj- The three spin-basis 
four- vectors s^, and form a right-handed system and provide, together with the 
momentum, an orthogonal basis system. They are chosen as: 

((p^o X Po, ) X p^H \ 
0, ,P , (B.4) 

''ix') = — (m,^Px^] ■ (B.6) 

rrij \ ^ \p^o\ ^ 



For a more detailed discussion see 14 
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B.2 Neutralino three-body decay — > Xk 

Here we give the analytical expressions for the different contributions to the decay 
density matrix for the three-body leptonic neutralino decay, where we sum over the 
spins of the final-state particles The contributions independent of the polarisa- 
tion of the neutralino x° 



= D{ZZ) + D{Z£l) + DiZin) + D{ldL) + DiidR) 



(B.7) 



are given by 
D{ZZ) 

D{Zh) 



9 



cos'^ 9w 



\/\{z)\\Li + 



X 



(B.8) 



9^ 



cos^ 6w 



U Re \ A(Z) fl^filA;{eL){20tj9i + Ot*9z) 



+ fl)fl,K{WOt;92 + Oyz) L (B.9) 



Di 



29' 



l/l/l/LP(|At(^~L)P^?i + |A„(£,)P^?2) 



+ Re{(/S)'(/L)%(^~L)A:(£^)}^3 



(B.IO) 



where A(Z') and At,u(fL) denote the propagators of the virtual particles in the direct 
channel and in both crossed channels (labelled t, m, cf. Fig. The quantities DiZi^ 
and D{ijii^) can be derived from Eqs. (|B.9|) and (|B.10|) by the substitutions 



The kinematic factors are 



^kj ^ ^kj ■> J Lj,k 



Rj,k ■ 



(B.ll) 



91 = {PxlPi-)iPx°Pe+) 

92 = {Px'iPi+)iPx''jP^-) 

93 = mjmk{pf~pf+) . 



(B.12) 
(B.13) 
(B.14) 



We can split the terms depending on the polarisation of the neutralino into T-even 
and T-odd contributions: 

S"z)(x?) = ^l^ix',) + ^^{X]) ■ (B.15) 
The T-even contributions depending on the polarisation of the decaying neutralino 



S^"(X^) = + T.''f{Z^L) + + + (B.16) 
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are 



T!'jf{ZZ) = 8 



cos'^ 9w 



|A(Z)p(i??-L^ 



^|^L,Re{A(Z) 
cos^ Ow ^ 



[(ReO,^/-(ImO,^/K + |0,^/K-^?2") , (B.17) 

'fyLkmL)i-2oj:^g'l + oj:;g^,) 

fL)fLK{h){20^,;9^, + O.Vs)] } , (B.18) 

l4fl/LniA«(4)P^?2"-|A,(4)lV] 

+ Re {{f%)\fl,f/^t{iL)K{h)gt}\ , (B.19) 



+ 



where the contributions lT^{Zi^ and {IbIr) are derived from Eqs. ( p.l8|) and 
( [B.19| ) by applying the substitutions given by Eq. ( [B.ll| ) and in addition (?! 2 3 
—gi2 3- The kinematic factors are 



9i = mj{p^o^pe-){pi+s'') , 

93 = mk[{p^ope+){pe~s'') - • 



(B.20) 
(B.21) 
(B.22) 



The T-odd contributions depending on the polarisation of the decaying neu- 



tralino 



ES^(X?) = SS"(ZZ) + + + ES'^(£i£i) + ES"(£^£^) (B.23) 



a,0. 



are 



cos^ 9w 

-A 



\^{Z)\\L\ - R\) 2Re {O^^) Im (0,^^.)< 



^1— L,Re{A(Z) 
cos"^ fc'vi/ 



+ /S/i,0|;.A:(4)]^4«}, 



= 2/ Re \ {fi*rVLrMh)Kih)9l 



l*\2( fl \2 , 



(B.24) 

(B.25) 
(B.26) 



where the contributions Yl'^'^{Ziji) and Yfjf^ {inlii) are derived from Eqs. ( [B.18| ) and 
( [B.19| ) by applying the substitutions given by Eq. ( [B.!!] ). The kinematic factor is 



gl = imke^,ypaS''''p'iop''g^Pi+ 



(B.27) 
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C. Kinematics 



C.l Phase space 

The complete cross section for the process can be decomposed into the production 
cross section and the branching ratios of the subsequent decays: 

datot = da{qg ^ q^ci) ■ dT{qL - qxl) ■ dT{xl - X?£+r ) , (C.l) 

where the factors E/mT come from the use of the narrow- width approximation for 
the propagators of the q^ and X2- This approximation is valid for (F/m)^ <C 1, which 
is satisfied for Fg^ ~ 4.7 GeV and for F^o ~ 3 x 10^^ GeV where the width is small 
because only the three-body decay is kinematically possible. 
We have: 

dT{qL ^ qxl) = -^P{xlq) d<^g , (C.2) 
dV{xl ^ xltn = ^D{xl) rf$^o , (C.3) 
where the phase-space factors in the laboratory system are given by: 

= 8{Lr ||p-,.,| cosS,. - 4- E,, - E,- coso:\ E,^'iE,^dn,Me- ,{C.5) 
where a is the opening angle between the two final state leptons. 
C.2 Integration limits 

When evaluating the phase-space integral at the parton level, kinematical limits 
need to be determined on the integration variables and these are listed below. The 
three-body decay phase space of the X2 has the following limits: 

ml - mj 

Ei- < — — ^ , C.6 

2(^x§ - \pS) 

2E^oEf- +m'? — mn 

cose,- < ' , ^ -. (C.7) 

2E,-|p,o| 



For a more detailed discussion see ||14| |. 
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